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We study the stability of patterns arising in rotating convection in weakly anisotropic systems 
using a modified Swift-Hohenberg equation. The anisotropy, either an endogenous characteristic of 
the system or induced by external forcing, can stabilize periodic rolls in the Kiippers-Lortz chaotic 
regime. For the particular case of rotating convection with time-modulated rotation where recently, 
in experiment, chiral patterns have been observed in otherwise Kiippers-Lortz-unstable regimes, we 
show how the underlying base-flow breaks the isotropy, thereby affecting the linear growth-rate of 
convection rolls in such a way as to stabilize spirals and targets. Throughout we compare analytical 
«j ' results to numerical simulations of the Swift-Hohenberg equation. 
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I. INTRODUCTION 



Pattern formation in thermal convection of a rotating fluid layer has been the subject of much experimental and 
theoretical work in recent years. The effect of the Coriolis force on the dynamics of thermal instabilities makes this 
system relevant for both astrophysical and geophysical fluid dynamics, while the appearance of spatio-temporally 
Oh \ chaotic dynamics near onset make it an attractive candidate for detailed analytical and numerical investigations of 
<—| ■ the origin and behavior of chaotic complex patterns. 

Kiippers and Lortz pj] determined that for rotation rates greater than a critical value Q > Q cr steady convective roll 
patterns are unstable to another set of rolls oriented at an angle j3 relative to the first. These results were confirmed 
and extended by Clever and Busse ^ who also determined the dependence of Vl cr and (3 on the Prandtl number of 
I . the fluid. In an infinite system, these dynamics are persistent due to isotropy. Busse and Heikes || used this fact, and 
the closeness of f3 to ^ to derive three coupled amplitude equations, in which rolls switch cyclicly as they approach a 
t-H ! heteroclinic orbit. In real systems, small amplitude noise perturbs this orbit, leading to nearly periodic switching of 
■ rolls. In sufficiently large systems the switching becomes incoherent in space and causes the development of patches 
\ of rolls with different orientations. The ensuing dynamics are chaotic, @], 0- 

In recent experiments on rotating convection ||, Thompson, Bajaj and Ahlers investigated the effect of a temporal 
modulation of the rotation rate on the Kiippers-Lortz (KL) state. They find that for sufficiently large modulation 
concentric roll patterns (targets) as well as multi-armed spirals can be stabilized and replace the chaotic KL state. 
Focusing on the target pattern, they find that the rolls in these patterns drift radially inward and they measure the 
dependence of the drift velocity on modulation amplitude and frequency, mean rotation rate, and heating. They point 
out, the modulation sets up an oscillatory azimuthal mean flow, which tends to align rolls along that direction. Since 
the alignment signles out a specific orientation, it breaks the isotropy of the system. Motivated by these findings we 
^ • therefore investigate here the effect of anisotropy on roll patterns in systems exhibiting KL chaos. 

Within the framework of a suitably extended SH-model we first study the stability of straight rolls in systems with 
broken chiral symmetry (modeling the Coriolis force due to rotation) and with weak anisotropy. We then use these 
analytical results to interpret simulation of this SH-model in a cylindrical geometry in which we obtain target and 
spiral patterns as seen in experiment. 



II. THE STABILITY OF ROLLS WITH ANISOTROPY 



We study the effect of weak anisotropy on the Kiippers-Lortz state in the following modified Swift-Hohenberg 
model, 

d t ih = fj,yj + a 2 (n- V) 2 ^ - (V 2 + 1) 2 V> - ip 3 + 7 fc • [V x [(W) 2 VV>]], (1) 

where n is a director indicating the preferred orientation, and a gives the strength of this anisotropy. We retain the up- 
down (Boussinesq) symmetry (ip — > —ip ) by including only odd terms in ip, and include a nonlinear gradient term that 
breaks the chiral symmetry. The rotation rate is therefore measured by 7. Similar models have been systematically 
derived from the fluid equations, with and without ^,|) mean flow effects, and have enjoyed widespread use, e.g. 
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JTo| O]. We mean (|T]) to be a model equation and are concerned with the qualitative effect of anisotropy on the 
Kiippers-Lortz instability. 

Focusing on the weakly nonlinear regime and assumimg the anisotropy to be weak, we take fi = e 2 fJ,2 and a = ea 2 
with e <C 1. To leading order in e the system is therefore isotropic. To study the effect of the anisotropy on the 
KL- instability we consider the weakly nonlinear competition of two sets of rolls with relative angle (3 with the ansatz, 



t/> = e(A(T)e t( - cos ^ x+sin< - e ^ + B{T)e l(cos{e+l3)x+sin(e+p)y) + c.c.) + h.o.t. 



(2) 



Thus here we do not analyze all side-band instabilities. To leading order the system is isotropic and 9 is a free 
parameter. The complex amplitudes A and B evolve on the slow timescale r = et. For concreteness we take n — i y . 
At order e 3 , a solvability condition yields 



d T A = n 2 A - a 2 sin 2 (9) A - 3| A\ 2 A - (6 + 4 7 sin f3 cos (3) \B\ 2 A, 
d T B = f i 2 B - a 2 sin 2 (6» + (3)B - 3\B\ 2 B - (6 - 4 7 sin/3cos P)\A\ 2 B. 



(3) 
(4) 



We examine the stability of rolls of orientation 9 with respect to a set of rolls oriented (3 to the first set of rolls. 
With a = (isotropic case), the absolute orientation of the rolls 9 is irrelevant, and we find they become first unstable 
to rolls with orientation 



KL 



45° 



for 



7 > 1KL 



Introducing a ^ leads to a dependence of both [3kl and jkl on the obsolute orientation of the rolls 
growth rates of the perturbations are given by 



a A = -2(/i 2 - a 2 sin 2 6), 
a B = A*2( — 1 + -7 sin /3 cos /3) — a(sin 2 ( 



(3+ [- 7 sin/3cos/3-2] sin 2 9). 
As can be seen from (^) the anisotropy has shifted the onset of rolls with orientation 9 to 

fj, 2 (0) = a\ sin 2 (9. 



(5) 

(6) 
The 

(7) 
(8) 

(9) 



Thus rolls with orientation 9 exist for ^2 > f L 2cr{9). For fixed fi this implies a neutral curve ct{9) as shown by the 
dashed line in Figure |l|. Rolls of orientation 9 first become unstable to rolls of different orientation at 



d/3 



0. 



(10) 
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FIG. 1. Linear stability diagram of rolls with orientation 6 in (|3|, |4|) with respect to rolls at a relative orientation of (3kl- 
Here fi — 0.2. Numerical results are given by the solid symbols: triangles for 7 = 3 and circles for 7 = 2. 
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Equation (|10|) must be solved numerically for the linear stability limits. Results are given in Figure |1| for various 
values of 7 (solid lines). To test these stability results we perform numerical simulations in which we perturb straight 
rolls of orientation 9 by small-amplitude rolls of orientation 9 + j3, where j3 is chosen as the angle corresponding to 
the maximal growth rate according to (|l0|). To verify that no additional instabilities are present, we also perturb the 
rolls with small-amplitude noise. As can be seen from the solid symbols in Figure [j], numerical simulations agree well 
with the weakly nonlinear analysis for rotation rates 7 that are not too far above "fKL{&) ~ 1.5, for which only weak 
anisotropy is needed for stability. For larger rotations rates 7, the weakly nonlinear theory overestimates the amount 
of anisotropy a needed to stabilize rolls. For a = 0(1) the anisotropy affects the linear growth rate of rolls already in 
(|l|) and will introduce a significant dependence of the critical wavenumber on the orientation 9. Numerical results for 
larger a reveal that large amplitude rolls (with 9 = 0) tend to grow and invade regions of rolls of other orientations 
front- wise. In fact, for a — > 00 only rolls with 9 = exist. 




FIG. 2. Stabilization of rolls in the regime of domain chaos arising from the Ktippers-Lortz instability, a) A typical patch-work 
pattern of domain chaos, where the angle between patches 13k l = 45, and 7 = 2.0, /1 = 0.2, a = 0.0). b) For the same values of 
the parameters with a 2 = 0.15, rolls are stabilized. 

Thus, weak anisotropy can stabilize periodic rolls arising in rotating convection in the Kiippers-Lortz unstable 
regime. Specifically, there is a finite band of angles 9 with respect to the anisotropic director n such that rolls with 
this angle are stable to homogenous perturbations of all possible orientations. 



III. MODULATED ROTATING CONVECTION: SPIRALS AND TARGETS 

We now turn to the specific problem of rotating convection with periodically modulated rotation. A thin layer of 
fluid of height d is heated from below and bounded above and below by a rigid plate, which is rotated with an angular 
velocity £1 = £1 (1 + Scosujt). For 5 = we recover the well-studied case of rotating convection §- §, Q^g. 
For i < 1 but nonzero, a nontrivial base flow is induced by the periodic motion of the rigid plates. This flow advects 
perturbations leading to thermal instabilities in such a way as to affect their growth rate. Indeed, far from the axis 
of rotation, the onset of the thermal instability is dependent on the orientation of periodic-roll perturbations with 
respect to the base flow in a manner analogous to the linear operator in (|l]), |L5|. Closer to the axis of rotation, the 
curvature of the base flow becomes significant and a straight-roll approximation is not a good one. 

The dynamics are given by the Boussinesq fluid equations in a frame rotating at the mean angular velocity Jl[- 
H with rigid boundary conditions, which require the azimuthal velocity component at the top and bottom plates to 
obey, 

u e = K(£O re ia;t ). (11) 

This condition induces an azimuthal shear flow, the strength of which grows with distance from the axis of rotation. 
If we assume our flow is restricted to a finite geometry, Coriolis forces acting on this flow can be balanced by the 
radial pressure gradient as with the centrifugal force. 

Far from the axis of rotation and, w.l.o.g., along the rc-axis, the non-dimensionalized flow takes the form, 

u e = Sm(prrx^ hkz -^ Z - 1) e iP "A , (12) 
\ 2 smli k I 



where the plates are located at z = 0, 1, k = + i) and 
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where the modulation frequency lu has been nondimensionalized with respect to the viscous diffusion time. Note that 
( |l2| ) satifies the equation of continuity, V • u = 0. 




FIG. 3. Representative patterns: a) KL-state with 5 = 0.0 b) 6-armed spiral for 5 = 0.0005 Equatii 
in this case, c) Target pattern for S = 0.001. For all three 7 = 2.0, fi — 0.2 and the system size L — 72. 




predicts six arms 



Thermal instabilities arising from the imposed temperature gradient are advected by (|T^) and hence paramctrically 
forced with frequency Proj. Rayleigh-Benard convection has been studied with time- modulated heating |l6|-|l9|], and 
time-modulated gravity (parametrically accelerated) [^0| . Here we note the following two facts which will allow us to 
deduce the effect of the base-flow ([jj]) on the instabilities to leading order. Firstly, near onset, the dynamics of the 
instability are slow compared with the period of the oscillatory shear-flow for any finite rotation rate. Secondly, the 
forcing is invariant under the transformation 6 — > —S,t — > t + p^j- Averaging over the fast oscillations with respect 
to the slow dynamics near onset, the base-flow affects the growth-rate of thermal instabilities through mean-squared 
contributions (proportional to <5 2 ). For the case (|l2|), the correction to the growth rate is of the form itg<9, 
consideration of the full, cylindrical geometry yields ^FUe(r) 2 dg. 

With the abovementioned assumptions, we study the dynamics of patterns in (fy) with, 



:2 , whereas 



= S 



V, 



(14) 



We note that the scaling of the anisotropy ( |l4| ) as linear in the distance from the axis of rotation, is only correct far 
from the axis itself. However, we retain this simplified form and hope to extract qualitatively correct results. In fact, 
simulations with other polynomial dependencies have revealed that only the monotonicity of the function is important 
in determining qualitative features of the patterns. 

Simulations reveal a wide variety of spiral patterns as well as targets. For small 6, where we expect the weakly- 
nonlinear theory for periodic rolls to be valid sufficiently far from the core of the spiral, we are able to predict the 
number of spiral-arms with reasonable accuracy. Such an analysis can be understood from Figure [l . For a fixed 
'rotation-rate' 7, the strength of anisotropy a increases with distance from the core of the spiral. There is thus a 
region in the vicinity of the core where no rolls are stable, and rolls of a given orientation 9* are selected at a distance 
r* as determined by condition (|l^). The number of resulting spiral-arms is then given by geometry as, 



N = qr* sin ( 



(15) 



where q is the wavenumber of the rolls (arms). 

Spirals or targets can be generated for the same parameter values given different initial conditions. In general an 
initial straight roll pattern will result in a target for sufficiently large 5 whereas disordered initial conditions gcncrically 
yield spirals, even for strong anisotropy. 

Interestingly, the orientation-selection mechanism given by ( p"5| ) predicts the possibility of a large region surrounding 
the core, within which no rolls are stable in the context of the weakly nonlinear theory. If the anisotropy is sufficiently 
weak, one should see a disordered region of domain chaos, bounded by a stable spiral, given a large enough system. 
Such a pattern is shown in Figure 0. 



IV. CONCLUSION 



Spirals and targets arising in Rayleigh-Benard convection have been the subject of much theory and experimental 
work, fip],fn]]2T|-|23[] . Target patterns in low-Prandtl number convection are a consequence of horizontal, thermal 
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gradients at the sidewalls of a cylindrical container, which tend to align rolls parallel to the walls [23|] . Even with 
sidewall forcing, the targets become unstable to straight rolls relatively close to threshold. In rotating convection, the 
target patterns arising from such sidewall forcing undergo a mean drift |ll]] due to the breaking of reflection symmetry 
by the applied rotation. However, in the case of rotating convection with a modulated rotation rate, the chiral patterns 
are not a consequence of the system geometry, but rather are induced by an isotropy-breaking shear flow, which acts 
azimuthally. We have shown that these patterns are stable in regimes where one would see spatio-temporal chaos in 
the absence of modulation. Our analysis indicates that the shear flow acts to stabilize rolls within a band of stable 
orientations w.r.t. the azimuthal flow itself. This leads naturally to a chiral pattern. 

The qualitative agreement between the types of patterns observed in experiment Q and those studied here make 
the selection mechanism described in section III plausible. Spirals and targets arise through the interaction of the 
destabilizing process responsible for the KL instability and the stabilizing effect of the azimuthal mean flow (MF). 
Quantitative comparison of the dependence of the pattern behavior on the reduced Rayleigh number, rotation rate, 
and amplitude and frequency of the modulation with experiment is, however, not possible within the framework of (Q). 
However, qualitatively the genericity of the appearance of spirals patterns under modulated rotation with disordered 
intial conditions (KL state) seems to hold both under experimental conditions and here. Targets, on the other hand, 
must be generated with care but then persist as stable patterns over a wide range of parameter values. 



FIG. 4. A stable spiral with a chaotic core. The system size is L = 144 with 7 = 2.0, fi = 0.2 and S = 0.0001. Note that 
equation ( |l5| ) predicts 14 arms. The core dynamics occur on a fast timescale w.r.t. the slow, solid-body rotation of the outer 
spiral. 

The target patterns observed in || travel inwards radially, collapsing periodically at the center. The origin of this 
drift has not been identified yet. One possibility is that the time-periodic component of the Coriolis force acting on the 
time-periodic azimuthal flow generates a radial flow with a steady component, which would advect the axisymmetric 
roll pattern [p4 25|. Another possibility is that the drift is due to a mismatch between the wavenumber selected by 



the umbilicus |26j and that selected by the container side- wall |27|,|2S|]. The competing selected wavenumbers set up 
a wavenumber gradient that induces a drift of the pattern |p9| , |30|| . It should be possible to distinguish between these 
two mechanisms by comparing the dynamics in systems of different aspect ratio. For larger systems, the wavenumber 
gradient induced by the incompatibility of the selected wavenumbers would be weakened, while the effective strength 
of the radial flow would naturally be stronger due to the greater Coriolis force at larger radii. Experiments in such 
larger systems would also be of interest in view of the prediction that in such systems the core of the spirals would 
exhibit chaotic dynamics of the Kiippers-Lortz type. 

In our simulations of (|l|) no radial drift of the concentric rolls was found. This is not unexpected since in the absence 
of the chiral-symmetry breaking term proportional to 7 eq.(]l|) is variational and persistent dynamics are ruled out. 
To obtain drift the variational character of the system has to be broken to a sufficient degree. This may require much 
larger rotation rates 7 or the introduction of additional non-variational terms. We have not pursued this further since 
simulations of extensions of (|l|) with possibly also modified boundary conditions would not allow any quantitative 
comparison with experiments and would therefore be of limited use in identifying the dominant mechanism responsible 
for the drift. 
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